Entanglement in Finitely Correlated Spin States 
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We derive bounds for the entanglement of a spin with an (adjacent and non-adjacent) interval 
of spins in an arbitrary pure finitely correlated state (FCS) on a chain of spins of any magnitude. 
Finitely correlated states are otherwise known as matrix product states or generalized valence-bond 
states. The bounds become exact in the limit of the entanglement of a single spin and the half-infinite 
chain to the right (or the left) of it. Our bounds provide a proof of the recent conjecture by Benatti, 
Hiesmayr, and Narnhofer that their necessary condition for non-vanishing entanglement in terms of 
a single spin and the "memory" of the FCS, is also sufficient jlj. Our result also generalizes the study 
of entanglement in the ground state of the AKLT model by Fan, Korepin, and Roychowdhury Q]. 
Our result permits one to calculate more efficiently, numerically and in some cases even analytically, 
the entanglement of arbitrary finitely correlated quantum spin chains. 
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Entanglement properties of quantum spin-chains have 
recently attracted attention from researchers in quantum 
information theory and condensed matter physics. From 
the perspective of quantum information theory, the dis- 
tribution of entanglement over long ranges via local op- 
erations on a spin-chain i3i,|4,,t5i| has obvious applications 
to teleportation-based models of quantum computation 
Moreover, it has recently been shown that en- 
tanglement in finitely correlated chains ^ can be used 
to achieve universal quantum computation |lO| and pro- 
vide a computational tool for adiabatic quantum compu- 
tation On the other hand, the scaling behavior of 
quantum correlations in infinite spin-chains is intimately 
related to their critical behavior (recent work has estab- 
lished a general mathematical framework for studying 
entanglement in infinite quantum spin-chains |12i|.) 

Motivated by the potential applications of distributed 
entanglement in finitely correlated chains, Benatti, et al. 
in 0, give a necessary condition for entanglement be- 
tween a spin and a subset of other spins; namely, that the 
entanglement between a spin and the "memory" of the 
finitely correlated state must be non-zero. They, further- 
more, conjecture that the same condition is sufficient, 
in the sense that it implies entanglement between a spin 
and a subset of other spins. We present here a proof 
of that conjecture by showing that the entanglement be- 
tween a spin and its neighbors converges exponentially 
fast (in the number of neighboring spins) to the entan- 
glement between a spin and the "memory" of the finitely 
correlated state. Moreover, we show that entanglement 
between distant spins vanishes exponentially fast in the 
length of their separation. 

Since finitely correlated states provide the exact 
ground states for generalized valence-bond solid mod- 
els our result generalizes the calculation of entan- 
glement for the AKLT model |Ti |. 

More importantly, our result implies a simple and com- 
putationally efficient way for detecting distributed entan- 
glement in finitely correlated states. Namely, the Posi- 



tive Partial Transpose (PPT or Peres-Horodecki) crite- 
rion can be applied to the state describing the 
interactions of a spin with the "memory" of the finitely 
correlated state, to detect entanglement between a spin 
and a subset of other spins. 



THE SETUP AND MAIN RESULT 

We will work with translation invariant pure finitely 
correlated states, FCS j^] on the infinite one-dimensional 
lattice. For each i € li, the spin at site i of the chain will 
be described by the algebra A dx d complex matrices. 
The observables of the spins in an interval, [TO,n], are 



given by the tensor product Au 



'l=„Mh- The 



algebra Az describing the infinite chain arises as a suit- 
able limit of the local tensor-product algebras A[-n,n] ■= 
®j=-n(-4)i- Any translation invariant state oj over Az is 
completely determined by a set of density matrices P[i,n], 
n > 1, which describe the state of n consecutive spins. 
In the case of a pure FCS, as was shown in these 
density matrices can be constructed as follows: 

The memory, B, of a FCS is represented by the algebra 
of 5 X 6 complex matrices. Let : A® B ^ B he & 
completely positive unital map of the form E(yl ® B) = 
V{A (g) B)V'', where V : (E) C'' ^ C^ is a linear map 
such that VV^ — i^B- We define the completely positive 
map E: B ^ B,hy E{B) = E(l^ (g> B). The condition 
on V implies that E is unital: E(lj3) = ts- For pure 
translation invariant FCS, one can assume that there is 
a unique, non-singular 6x6 density matrix, p, such that 
TrpE(B) = TrpS, for aU S e 6. 

We introduce the density matrix pAt^B associated with 
the state encoding the interaction between the spin at site 
1 and the "memory" of the FCS: 



Tr^^B [pA®B A®B^= Tre [pE{A® B) 
Using the cyclicity of trace we also have PA®B = V'^ pV . 
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We are now ready to define the density matrices P[i.n]- 

where F„ = (1^ ® F)(l^»2 (g) • • ■ (g) F). 

An important property, intimately related to the ex- 
ponential decay of correlations in a pure FCS is that the 
peripheral spectrum of E is trivial; that is, Ig is the only 
eigenvector of E with eigenvalue of modulus 1 0] . This 
implies that the iterates of E converge exponentially fast 
to E°° given by t°°{B) = lim„^ooE"(B) = Tr(pB)lg. 
More precisely, for any A such that |Ai| < A < 1, for 
all eigenvalues Xi of E different from 1, there exists a 
constant c such that for all n > 1: 

||E" -E°°|| < cA", (1) 

where the norm is the oo-norm on B considered as a 
Banach space with the 1-norm. 

Our object of study is the entanglement of formation, 
EoF The EoF is defined for states of composite 

systems with a tensor product algebra of observables 

Definition 1 (Entanglement of Formation). The en- 
tanglement of formation of a bipartite state over Xi 
with associated density matrix ai2 is given by: 

i 

where S{p) = — Trplogp is the von Neumann entropy 
and the infimum of the average entropy is taken over all 
convex decompositions oyi = ^iPi '^12 ''^^'^ pure states. 

Whenever Xi is finite dimensional, as will be the case 
for us, the infimum can be replaced by a minimum in 
the above definition, i.e., there is an optimal decompo- 
sition, {pi, a\2}, where the infimum is attained (see [l9j 
for details). We call an ensemble for the den- 

sity matrix a whenever the latter can be decomposed as 
(7 = '^^\4'i){4'i\- There are an infinite number of en- 
sembles corresponding to a given density matrix. The 
following lemma provides us with a complete classifica- 
tion: 

Lemma 2 (Isometric Freedom in Ensembles, [20|). 

Let {\ei)}f^i be the ensemble corresponding to the eigen- 
decomposition of the density matrix a, where d = 
rank(a). Then, *s ensemble for a if and 

only if there exists an isometry [/ : t-> C™ such that 

d 

IV'i) = XI l*^-?'^' 1 < « ^ 

The above lemma implies that any two ensembles for 
the same density matrix, )}*=!, and are 
similarly related via a partial isometry W : C*^i C*^" . 
Our main result is the following theorem: 



Theorem 3. For any pure translation invariant FCS we 
have 

< Eij^^s]iPA(g,B) - E[_A,J^(s,^-l]{p^l,n]) < e(n), (2) 
where e{n) decays exponentially fast in n. 

PROOF OF THE THEOREM 

The lower bound is proven in 01 . For the sake of com- 
pleteness, we include here the following proof. 

The definition of implies that every decomposi- 

tion of pAigiB into pure states induces a decomposition 
of P[i,„]. Moreover, the restrictions to the spin at site 1 
of the i-th state in the corresponding decompositions of 
PA®B and P[i,Ti] a-re equal. To see this, note that since the 
operators Vn leave the first spin invariant, the cyclicity 
of the trace implies 

Tr^»"-i(P[i,„]) = Tr^»"-i®i3(Klp!4<»B^") ^ TrB(p^g,g), 

where we have used VnV^ — tj( ® 1b- It follows that 
for each decomposition of pj^tsB there is a correspond- 
ing decomposition of P[i^n\ with equal average entropy. 
Since the average entropy of p[i^n\ is minimized over a 
(possibly) larger set of decompositions, the lower bound 
follows. 

We now focus on the upper bound. We start with the 
following decompositions of pa®b and p[i.„] into (unnor- 
malized) pure states: 

6 

PA^B = (3) 
1=1 
b 

= E Gl,y^\x^){X^\VG^,3, (4) 

where {Ixi)},^^! is the eigen-ensemble of p and Gn,j = 
® IXi)/IIXjll)- The term in parenthesis in the 
expression for G„j comes from the Kraus operators in 
the decomposition of the completely positive map Trg. 

By the observation following Lemma 13 we have that 
the (unnormalized) states |<1>") in the optimal decom- 
position of are given by: 

b 

I^D- E 1<1<L, (5) 

,2 J 

for some partial isometry U : <C i-^ C , whose depen- 
dence on n we suppress. Moreover, it is easy to check 
that: 

b 
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is an ensemble for 

To calculate the EoF we need the restrictions of 
and{\^i){^i\}to A: 

= Tr^«.-i(|<fr)(*ri), i'l = (7) 

Define the density matrices — 0" /a" and '4'i = ij^i /Pi , 

where = ||^rlli = M4>?), A = WMi = Tr(?A/). 
From the definition of the EoF and the optimality of 
we get: 

L 

E[A,B(g>8]{PAm) ~ E[y^ j^i)n.-i]{p[i n]) < ^e/(n), (8) 

1=1 

where ei{n) = PiS{^i) - afSic^]"). 

It remains to show that X^zLi is exponentially 

small. We estimate each term in the sum as: 

IqWI <AI^(^0- W)l + IA-«riiogd, (9) 

since rank(0") < d. 

To bound \S{ipi) — S{(j)")\ we use Fannes' inequality for 
the continuity of the von Neumann entropy [2 If : 

\S{i;i) - SW)\ < {\ogd + 2)Ui - 0rili +'7(11^; - 

where ri(x) — — a:logx and log is the natural logarithm. 
By the triangle inequality we have: 

i/3,-ari = iii^dii-ii0riiii<ii^;-^riii- (n) 

Another application of the triangle inequality gives: 

\\Mi~aY<l>ni + \\K-f3i)'P?\\i 



I - Vi 111 



< 



which simplifies, with the use of to the following 

inequality: 



7 ^ (Pi 111 



< 2- 



Pi 



(12) 



Combining equations (|^- (|12|l and setting 

rr^Ui-m\ilPu (13) 

we get the following bound for e/(n): 

\ei{n)\ < A[(logd3 + 4)rr + r;(2rr)]. (14) 

where we have assumed that 2r" < 1/e, to assure ri{x) 
is increasing. 

To complete the proof, we show that r" is exponen- 
tially small for large n. Since each Gnj leaves the spin 
at site 1 invariant, the cyclicity of the trace yields: 



ButG„,,<,, = l^®IE"-i(|x,)(x,'|)/(||x,|l||x,'||). Sub- 
stituting E'^^ for E"^^ we get: 

b 

i,*' Jj' = l 

where X,^,, = V^Xt){xe\V and Y^^j, ^ 1a ^ [E""' - 

E°°](|X.)(X,'I)/(IIX,IIIIX,'II)- 

Like all trace preserving quantum operations, the par- 
tial trace is contractive with respect to the 1-norm. 
Hence, an application of the triangle inequality for the 
1-norm gives: 

b 

||^;-0rili< E |C^;V/)llf^^fo)lllAHIirw'lli 

hi',j-J' 

It is not hard to see that 

II^M'lli^llx.llllxdl, r,,,'||i<l|]E("-i)-E-|| 
and hence 

11^; - ^riii < (E|E I iixdi)' iiE^"^'^ - E°°ii 
1=1 j=i 

Since |C^;,(y) P ||x»lP = Ph two applications of 

Cauchy-Schwarz give: 

ii^i-^nii <&'aiie(""'^-e°°ii. 



(15) 



Finally, combining ^ with equation l(T^ be- 

comes: 

T," < ciA", ci = cb^/X. (16) 

To conclude the proof, we note that since the bound for 
t" is independent of I, summing over I in equation l|14|) 
yields: 

L 

Ekn < (logd=^+4)ciA"+r7(2ciA"). 

1=1 

It is clear that for A' > A there exists a constant C2 
such that 

??(2ciA") < C2(A')". 

The only condition on n was imposed in equation p4|) 
were we assumed that 2t" < i. Using equation 11()|) we 
see that there is an uq such that the above condition is 
satisfied for all n > hq. The previous observations imply 
that for all A' with A < A' < 1, there is a constant C3 
such that: 



E Ki^'J')ULi^n^^MV^\X^){X^'\VG,,^,Gl^,), 



e(,^)=C3(A')">El^ri, for all 



1=1 
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Finally, equation (O implies that: 

EiA,Bm]{PA<g)B) - -E[^,^®"-i](p[i,„]) < e{n), 

and this completes the proof of the theorem. 

A natural question to ask at this point is the following: 
How does the entanglement between the spin at site 1 and 
spins at sites [p,n], p > 2 behave as p becomes large ? 
Since the state Pi,[p,n] factorizes into pi ® p[p,„] as p ^ 
00 [i^l, we expect that the bulk of the entanglement is 
concentrated near site 1 . The following theorem confirms 
this: 

Theorem 4. For any pure translation invariant FCS and 
n>p>2, the following bound holds: 



where e(p) decays exponentially fast in p. 



(17) 



Sketch of the proof: The main observation is that the 
trace distance between the states Pi.[p,n] and pi ® P[p,n] 
vanishes exponentially fast with p. This is a consequence 
of the exponential rate of convergence described in equa- 
tion (Q. Since E]^a,A»^-p+^{Pi ® P[p,n]) = 0, a straight- 
forward application of Nielsen's inequality for the conti- 
nuity of the EoF [2^ yields the desired result. □ 



DISCUSSION 

Having established such a strong connection between 
the states P[i^n\ and pa<^B, one can apply various entan- 
glement criteria on Pai^B to deduce entanglement prop- 
erties of the spin chain. To start with, we note that 
for qubit chains with 2-dimensional memory algebra i3, 
the entanglement of P[i^n] can be computed analytically 
(in the limit) by evaluating the concurrence ji^l oi pA®B- 
For higher dimensions one can apply the PPT criterion to 
PA®B to detect distributed entanglement in the fini tely 
correlated state. Specifically, the main theorem in |25| 
implies that there can be no PPT bound entanglement 
in pA<sB since rank(p_4^g) = b < maxjd, 6}. Hence, if 
the partial transpose of pa(sb is positive, then pAt^B is 
separable. On the other hand, if the partial transpose 
of PA(SB is negative, then for n large enough P[i^n] be- 
comes entangled. The amount of maximum entangle- 
ment in P[i^n] depends on the amount of entanglement 
found in PA®B- From this point of view, it would be very 
interesting to look at FCS that maximize entanglement 
of Pai^B- Moreover, understanding how entanglement of 
PA®B varies with different CP maps E could lead to a bet- 
ter understanding of how phase transitions occur when 
we vary the parameters in the underlying hamiltonian of 
the system. 

To conclude, we note that the conjecture of Benatti, et 
al. 0, follows as a corollary of Theorem|2| In particular. 



our result implies that a spin at site 1 of the chain is 
entangled with spins at sites [2,n] (for n large enough) 
if and only if PA®B is entangled. Moreover, the entan- 
glement of P[i,Ti] approaches the entanglement of pa®b 
exponentially fast. 
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